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We study the shape and texture of finite domains comprising chiral/achiral molecules carrying 
a tilt field embedded in a 2-dimensional surface. Using a combination of simulations and exact 
variational calculations, we determine the equilibrium (mean field) phase diagram in the achiral and 
chiral cases. We find a variety of novel shapes and textures including a spontaneously broken chiral 
texture when the molecules are achiral. We show that chiral tilt-domains nucleating in a region of 
two-phase coexistence, repel each other, thereby preventing coalescence and further growth. Our 
work has implications for tilt texture domains in phospholipid monolayers and giant unilamellar 
vesicles (GUVs), nucleating domains of Sm-C* in Sm-A films, and chiral emulsions in Sm-A films. 



INTRODUCTION 



Shapes and textures of finite-size domains embedded in two or three dimensional substrates are governed by 
an interplay between elasticity, defects and finite size. Such self-assembled patterning gives rise to unique 
rheological and optical properties with potential engineering applications as in the production of templates, 
micro-lithography, liquid-crystal displays and optical switches [lj. The purpose of this paper is to provide 
a detailed analytical and numerical study of the shapes and textures of finite-size domains in 2-dimensions 
when the constituent molecules (or molecular aggregates) are both chiral and possess a tilt 0, Q. This 
has implications for a wide variety of soft materials such as (i) domains in freely suspended phospholipid 



monolayers in an air-water interface 0, [H, 0, B, ( ii) domains in tense lipid bilayer vesicles [1, [§, E3| or 
freely suspended bilayers in an aqueous solvent (iii) nucleation of Sm-C* domains in freely suspended 
Sm-A films at the A-C* transition [l2j], (iv) droplets of a suspension of ferromagnetic particles in a liquid 
film [131 ] and (v) emulsion of chiral nematogens in a liquid film. A variety of techniques such as fluorescence 
confocal polarization microscopy, atomic force microscopy, Brewster-angle microscopy and Bitter microscopy 
(in the case of ferrofluids) have been used to detect domain shapes and textures on 2-dimensional substrates. 
Moreover the tilt phase of chiral molecules with an electric dipole moment is ferroelectric [lil . [lj| , allowing 
for the possibility that external electric fields may be used to bias the textures 

Over the years there has been a substantial amount of experimental and theoretical work in this area. 
However, most theoretical studies have been restricted to the case when the constituent molecules are achiral. 
In this paper we systematically study the effects of chirality on the shape and texture of tilt domains. The 
competition between bulk orientational order, chirality and anisotropic line tension from the boundary gives 
rise to a variety of shapes and textures. Even in the achiral limit we point out several new phases missing in 
earlier analyses, such as the annular phases, which we establish using an 'exact' variational approach (whose 
meaning will be made clear later on). Unlike previous studies, our variational calculation allows the topology 
of the domain to vary, i.e., the domain need not be simply connected. The method of analysis followed in 
this paper takes off from and extends the elegant treatment of Ref . [jlj . 

In the next section, we discuss the various contributions to the domain energy containing chiral molecules 
with an orientational field. We will arrive at an effective continuum hamiltonian (Sect. 2.1) describing the 
shape and texture of this finite, two-dimensional domain. We will next discuss (Sect. 3) the optimum shape 
and orientational ordering within the finite patch of fixed area. The main results of this paper are — (i) We 
obtain an 'exact' phase diagram in the achiral limit when the boundary of the domain is circular (Sect. 3.1), 
with a variety of optimal textures such as virtual boojum, defect, achiral annular, and a spontaneously broken 
chiral annular phase, (ii) We study smooth perturbations of the boundary, determine the regime of instability 
of the circular domain (Sect. 3.2), and the nature of the resulting equilibrium phase. We find that the annular 
phases are no longer minimum energy configurations. Typical domain shapes are thin and elongated, best 



described as being rectangular with circular caps, (iii) Turning on bulk chirality gives rise to a spiral defec? 
phase and a novel chiral tweed phase (Sect. 3.3). (iv) Chirality may even induce a large enough domain to 
split into multiple domains; we determine the conditions under which such multiple domains obtain (Sect. 4). 
This suggests that tilt-domains, nucleated following a quench across a phase-boundary in a multicomponent 
system, would repel each other, preventing coalescence and subsequent growth. Finally we discuss how an 
electric field might couple to the tilt when the molecules are chiral; this may be used to switch between 
degenerate textures (Sect. 6). In the following sections we will see how these results are obtained. 



II. CONSTRUCTION OF AN EFFECTIVE HAMILTONIAN 



The question we address is : what is the equilibrium texture of a collection of chiral molecules described by 
an orientational field, uniformly spread over a finite domain of fixed area A embedded in a 2-dimensional 
(2d) flat substrate ? We show that the interplay between orientational ordering and domain shape gives rise 
to a variety of novel phases and shape transitions. We will first discuss the nature of the order parameter 
and the dominant contributions to the energy of the domain. 

The orientation of a rigid molecule may be described by a polar vector which takes values in S 2 (Heisen- 
berg spin) or a director which takes values in the projective space RP 2 = S 2 /Z2 |l8|. However energetic 
considerations, e.g., in the case of membranes, a combination of van der Waals and hydrophobic shielding, 
may constrain (a) the centre of mass of the molecules to lie on a 2d, flat surface and (b) the projection 
of the long-axis of the molecule onto the 2d plane to have a fixed magnitude (or small deviations from a 
fixed value). Thus owing to strong uniaxial anisotropy, the low energy sector may be described by a 2d 
polar vector m (since the 2d surface carries a unique local outward normal) which takes values in S 1 , an 
XY spin. The centre-of-mass density p(x, y) of the molecules is assumed to be uniformly smeared over the 
patch of area A. In addition, individual molecules may have a permanent dipole or quadrupole moment — 
this is especially true in the case of zwitterionic lipids (electric dipole moment) or ferromagnetic particles 
(magnetic dipole moment). The molecules outside the finite patch are assumed to be either in the isotropic 
or the Sm-A (i.e., no tilt) or the liquid-disordered (Ld) phase. 

The tilt molecules interact with each other, and with the molecules outside the patch, both sterically (purely 
repulsive) and via short range (e.g., van der Waals) attractive interactions. Both these effects contribute to 
chiral interactions; the former via the Straley picture of interlocking screws [19], the latter via a generalisation 
of the Van der Waals dispersion to chiral molecules 20]. In the continuum limit, these short-range interactions 
can be written as the usual Frank energy, modified to include the effects of chirality. In addition to these 
short-range interactions there could be long-range dipole-dipole (or higher multipole) interactions between 
the tilt molecules carrying a permanent dipole moment, which may also contribute to chiral interactions. 

In this paper, we ignore the contribution coming from dipolar interactions, which would be justified if they 
are demonstrably smaller than the Frank terms (see Sect. 5). For instance, in the case of phospholipid 
domains, the Frank elastic constants arc of order 10fcsT room , while the dipolar energy per phospholipid is 
about 10 -22 J, and hence may be ignored. In situations where dipolar energy scales are comparable, it is 
likely to have a major effect on domain shapes and textures, such as, for instance, the formation of string-like 
aggregates [2l|. 

In our present treatment we assume the substrate to be flat over the scale of the patch; thermal undulations 
of the substrate (liquid film or membrane) controlled by surface tension or bending rigidity are assumed 
to be negligible over this scale. Elsewhere we have discussed the effect of chiral tilt textures on a flexible 
membrane, leading to a novel chirality induced budding and tubulation (23l . l24j . 



A. 



Tilt texture hamiltonian 
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Let us assume that the molecules, now described by a 2d orientational field m, are smeared with a fixed, 
uniform density over a patch (domain) of area A embedded in a 2d surface, the xy plane. The rest of the 
surface surrounding the patch is a structureless fluid consisting of molecules different from the constituents 
of the patch (or else in a different phase) or a solid substrate. We assume that the patch constituents have 
come together as a result of micro phase segregation. The domain energy can be written in terms of bulk 
distortions of m and interfacial (perimeter) distortions, including an orientation- interface coupling. 

In our low energy description, the magnitude of m is held fixed and normalised to unity (hard spin model). 
The resulting domain energy only contains phase distortions of m = (cos 0, sin 0). Such a phase- only theory 
is formally ultraviolet divergent, and needs to be made finite by a microscopic cut-off length of molecular 
dimensions. The regime of validity of such phase-only theories is over scales larger than this cut-off length. At 
these scales, the hard-spin model is equivalent to a soft-spin model, where the magnitude of m is constrained 
by means of a potential ^(m) = u (m • m — l) 2 to deviate only slightly from its preferred value, for u large 
and postive. 

The form of the bulk distortion energy of m can be constructed from symmetry arguments. Since tilt is 
a spatial vector, tilt distortions may be described by a hamiltonian invariant under 2d spatial rotations 
0(2). To lowest order, this leads to the usual Frank-energy in 2d, with two independent contributions - 
a splay distortion (V • m) 2 and a bend distortion (V x m) 2 , where V = (d x ,d y ). Recall however that the 
molecules are chiral, and so on general symmetry grounds, we should allow for terms which are invariant 
under improper rotations alone; thus to lowest order this gives rise to a term of the form (V • m) (V x m). 
Note that for a 2d vector field m, V x m is a pseudo-scalar, not a vector. Viewed as a vector in R 3 , the curlm 
is a pseudo-vector pointing along z (we will consistently represent 2d vectors by boldface and 3d vectors by 
an overhead arrow). 

The coupling between the tilt and interface gives rise to boundary terms — if n is the local unit normal to the 
boundary aiming into the domain, then there are two possible boundary terms to linear order in m. These 
are (i) an isotropic line tension proportional to the perimeter and (ii) anchoring terms of the form (m • n) 
and (m x n), which contribute to an effective anisotropic line tension. One source for this boundary term 
is that the interface between the tilt and non-tilt components (in the case of chiral emulsions this coincides 
with the tilt-solvent interface) prefers a specific alignment of the tilt with the local normal. Alternatively 
such boundary contributions may arise from the presence of a bulk spontaneous bend V x m and splay V • m, 
since 



The spontaneous bend is a chiral contribution; spontaneous splay could arise from specific substrate-tilt 
coupling or from steric forces arising from large head-to-tail ratio of the constituent molecules. 

Thus our low energy effective hamiltonian for the chiral tilt domain in the strong segregation limit, is given, 
to quadratic order in the fields, by E = Eb + Ec, where the bulk energy 
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and the interfacial energy, 
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When the coefficients k\ and ki are comparable, the Frank terms resist any deviation of m from uniformity 
and hence the m field would everywhere point in a specific direction. The anisotropic boundary terms 
however prefer to align m at the boundary along (or orthogonal to) the local normal. This competition gives 
rise to a variety of nontrivial textures. The boundary terms may be interpreted as providing an effective 
anisotropic line tension of the form cto + <j\ (m • n) + 02 (m x n). A negative value of this effective line tension 
would lead to an instability of the circular domain boundary [llj • Nontrivial textures also obtain when the 
Frank coeffcients ki and ki are appreciably different [2^]. We will see that chirality introduces a new twist 
to the texture phase diagram [23, [24j . 

In this strong segregation limit, molecules are not allowed to exchange across the domain. This imposes a 
constraint on the allowable configurations explored by the molecules. Thus for fixed "bath" conditions (state 
of the molecules outside the patch), the configuration of molecules inside the patch can reach equilibrium. 
The area of the patch and the density of molecules in the patch are fixed over this time scale. In the case of 
nucleation of the chiral tilt domains (e.g., Sm-C* domains in a freely suspended Sm-A film or tilt domains 
of chiral lipids in a mono/bilayer), slow domain growth justifies the use of quasi-equilibrium ideas. For a 
dilute emulsion, the area of the domain is fixed at time scales smaller than domain coalescence times; thus 
the shape and texture of the domain will be such as to minimise the free-energy subject to the constraint of 
constant area. 

Our aim then is to find that optimal conformation, in general a difficult calculational task. We will see that 
given the energy scales in the problem (Sect. 5), we may ignore the effects of thermal fluctuations on the shape 
and texture of the domain. Thus, the Frank constants ki, ki are ~ 10/cb T room while the interfacial energy 
of a domain of radius 10 fim is an order larger. Hence a mean field minimisation of the above free-energy 
functional will suffice. 



III. MEAN FIELD PHASE DIAGRAM 



We will perform this constraint minimisation of the free energy using a variational approach, where we vary 
both the texture and the domain shape while keeping the area constant. Our constrained variational ansatz 
also accounts for the possibility that the domain may not be simply connected, however we will not allow 
the domain to break up (domain splitting will be treated in Sect. 4). Our variational guesses are supported 
by computer simulations, i.e., monte carlo simulations with simulated annealing to avoid getting stuck in 
local minima. It is in this sense that our variational calculation for the mean field phase diagram is 'exact' 
(made more explicit in the next subsection). We have also explicitly checked that the effect of including 
higher order (symmetry allowed) terms in the hamiltonian is small and does not affect the phase diagram. 

We henceforth set ki = 1 as our unit of energy and R = 1, associated with the domain size R — \J A/n, 
as our unit of length. Some experimental situations are more conveniently analysed by fixing <jq — 1 as the 
unit of length and determining the phase structure in the (R, a) plane. We discuss both these ensembles. 



A. Achiral domain 



We first turn off the effects of bulk chirality (k c = 0), but retain the anisotropic boundary term. In most of 
our work we take ki = k\ = K (the one-coupling-constant approximation); towards the end we comment on 
the distinct features when k\ ^ ki. 



FIG. 1: Effect of anisotropic line tension on the texture of an achiral circular domain : (A) a boojum of charge 2, 
(B) texture of a domain of radius R whose centre is at a distance a from the core of the boojum. The arrows give 
the local direction of m. 



In terms of the phase angle 4>, where m = (cos 0, sin 0), the free-energy functional simplifies to 

E = — J (V0) 2 + j> dl (a + a x cosO - 0) + cr 2 sin(i/> - 0)) (5) 

where we have written the normal n = (cos tp, sin ip) in polar coordinates. Variation of the above hamiltonian 
gives the bulk Euler-Lagrange equation, V 2 = (Laplace) and nontrivial boundary equations which have 
to be simultaneously satisfied. Instead of asking for solutions of the combined bulk+boundary equations, we 
determine the global minimum free-energy configuration variationally — this has the advantage of obtaining 
boundary minima, minimum energy configurations that may not be solutions of the Euler-Lagrange equations. 



Without loss of generality we may work with a simplified free-energy functional — as argued by previous 
authors [T3, G3]) a global rotation of m by (p' = <p — tan -1 ((72/01) — ir/2 without rotating the boundary 
contour transforms the free-energy functional to a new functional 

E = y J (V0) 2 + j dl {ct + a cos(0 - i[>)) (6) 

where a = \/ o~\ + o\ is the single redefined anisotropic boundary coefficient. Configurations which minimise 
the free energy ([5]) with parameters (K, o~q , o\ , o~i ) can be obtained from configurations which minimise the 
free energy ^ with parameters (K,a ,a) via the global transformation defined above. Thus in this equal- 
constants, achiral regime, it suffices to work out the phase diagram in the (cr 0l a) plane using the free energy 
(frj|), the phase diagram for arbitrary values of <7i, o~i can then be reconstructed. This simplification does not 
occur when the bulk chiral term is present. Note that by making a global rotation, it is possible to entirely 
climate the contribution from the boundary chiral term V x m. This implies that the domain shapes will 
always be achiral fl7j | . 



Now set R = 1, as our unit of length. We will first discuss the limit of large oq, when the domains are forced 
to be circular (we will discuss the limit of stability of a circular domain in the next section) . The extreme 
limit, (To — > oo, might be experimentally arranged by smearing a thin film of Sm-C* over a hole made in 
a solid substrate. Clearly in the absence of any anisotropic line tension, a = 0, the lines of m describing 
the optimal texture are a (parallel) set of straight lines pointing in any direction. This configuration is 
invariant under arbitrary translations and rotations of the boundary circle with respect to the lines of m. 
An infinitesimal (positive) a forces the lines of m to curve slightly (at the cost of Frank energy) to meet 
the circular domain boundary at a desired angle (Fig. 1). From the axial symmetry of the solution, it is 
clear that this configuration is associated with a +2 virtual defect, situated outside the domain, called a 
virtual boojum [12j. Note that this configuration is infinitely degenerate. To calculate the energy we need to 
parametrize the virtual boojum configuration, where the centre of the circular domain is a distance, a, away 
from the core of a virtual defect of strength N. Placing the origin of polar coordinates (r, 9) at the core of 
the boojum, we can describe the texture m = (cos 4>, sin (f>) by the equation 



4> = NO + c\r sin 9 + c 2 r 2 sin 29 + c 3 r 3 sin 39 + 



(7) 
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FIG. 2: Variation of a, the distance of the virtual boojum from the centre of the domain, as function of the anisotropic 
line tension a : note that the penetration of the boojum core into the domain is pre-empted by the annular phase. 
We have taken ctq = 1 for convenience. 



The first term is the singular part representing the boojum, subsequent terms represent the most general 
smooth solution of Laplace's equation with the desired symmetry. The energy E is minimised with respect 
to the variational parameters (N, a, {c„}) for fixed (o"o,er). Consistent with our symmetry arguments, we 
explicitly find that the optimal value of N is exactly 2 and independent of a. In principle ofcourse we can 
take as many variational parameters {c„} as we desire, however in drawing the phase diagram we have 
parametrised the texture by a and c\ alone; inclusion of higher order c n lowers the energy by 1% at most. 
The form of the variational energy for (j) = 28 is given by 



E(a) = -2irK\n 



. i * 

a 



27T h 27TC70.R 



(8) 



which is to be minimised with respect to a. The minimisation with respect to the whole set of parameters 
(a, {c„}), however, is best carried out numerically. Upon variation, we find that an increase in the anisotropic 
line tension a leads to an increase in ci, pushing the core of the virtual defect towards the domain centre. 
Experimental characterisation of this texture in the context of Langmuir monolayers at the air- water interface 
0, 0] agrees with our analysis above. This movement of the core continues till a = 1.7, beyond which there 
is a discontinuous transition to the annular phase (Fig. 2). Note that the boojum is always virtual, the core 
never penetrates the domain, since it is pre-empted by the annular phase. Reference [l7j misses this feature, 
since their variational shapes do not include such multiply connected topologies. 



The Achiral Annular phase can be parametrized by an inner and outer radii r\ and r-i respectively. At the 
outer rim of the annulus, m is directed radially outward, while at the inner boundary, m is inclined at an 
angle a to the local normal n. With the origin of polar coordinates at the centre of the annulus, the texture 
may be described by 



r-2 



r2 - n 



(9) 



where r\,a are variational parameters (r2 may be obtained from the constant area constraint). With this 
parametrisation, the energy of this configuration is given by, 



E(r\, a) = ttK 



\a 2 


( T2 


+ n\ 


+ ln^" 
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27Tcr (r2 — ri cos a) + 27rcro (n + r 2) 



(10) 



where r\ — r 



R 2 . The inner radius r% monotonically decreases as we increase a (Fig. 3) such that m 
always points radially outward (cv = 0) at every point of the domain. With increasing a, we cross a first 
order phase boundary (Fig. 4) into a new annular phase with a chiral texture. 



The Chiral Annular phase obtains when the angle a jumps to a value greater than 7r/2 at both the inner 
and outer boundaries. The bulk texture is chiral and is a true symmetry broken configuration (recall that 
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FIG. 3: The inner radius n of the annular domain decreases (increases) with a in the achiral (chiral) phases (repre- 
sented by circles (squares)). We have taken ao — 1 for convenience. 



there is no bulk chiral interaction) with the pseudoscalar order parameter 



This chiral phase is doubly degenerate and can spontaneously acquire either sign; the anisotropic tension a 
behaves as a surface field conjugate to this Ising-like chiral order parameter. 

Note that the 'core' of these annuli must consist of molecules without tilt; our assumption in minimising the 
energy has been that the two species of molecules are free to diffuse across the domain of fixed area A. If 
for some reason the diffusion of the non-tilt molecules from outside the patch to the 'core' of the annulus is 
hindered, then the appearance of such a phase would be kinetically blocked. In this case, the virtual boojum 
will penetrate the domain and would lead ultimately to the Hedgehog configuration (as in [13 )• 

The Hedgehog phase has a defect at the domain centre and a texture described by (f> — 9 (when the origin 
is at the domain centre). The energy includes a core contribution e c of the defect of microscopic size r c 



The phase diagram in the (R, a) plane is shown in Fig. 4. The transitions between these phases are discon- 
tinuous, as indicated by the change in slope of the energy branches as a function of a. As discussed earlier, 
in some experimental situations it is more useful to study the phase diagram in the (R, a) plane; thus we set 
(To = 1 to obtain a unit of length. Assuming that the domain is still circular, the phase diagram shows the 
discontinuous transitions which weaken (i.e., the jump in the appropriate order parameter decreases) as the 
domain size R shrinks. 

This analysis suffices to reconstruct the entire phase diagram in the {a 1^2) plane by global rotation of 
optimal textures obtained above. For instance, the virtual boojum texture reported in [171 ] can be obtained 
from a global rotation of the texture in Fig. 1. Likewise the Vortex phase can be obtained by globally rotating 
the hedgehog texture by an angle tt/2. The vortex texture spontaneously breaks chiral symmetry. Other 
nontrivial phases may be obtained by global rotations of the texture in the annular phases. 

Finally we mention in passing, the case when k\ ^ k^, which has been the subject of study in [22|. The 
difference in the splay and bend energies implies that the hedgehog and vortex defects are no longer degen- 
erate. If AK — (ki — k2)/ki, then as AK goes from large negative to large positive values, the texture goes 
from a hedgehog to a vortex via a spiral defect phase. 




(11) 



(r c < R), 



Ehedge = kK In 



2irR (a — ctq) + e c 



(12) 
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FIG. 4: Phases of the tilt texture domain with circular periphery : (1) Virtual Boojum (2) Achiral Annular (3) Chiral 
Annular (4) Hedgehog with e c = 0. Arrows indicate direction of m field. 



For smaller values of <7q/ Vi, the domain is no longer circular; as shown by [17ll25|. deviations from circularity 
arise when the effective line tension cto + 0i(m • n) + 02 (m x n) < 0. We summarize their results, based on 
linear stability analysis. 

Parametrize the boundary by smooth deformations of a circle, r{6) where we have temporarily shifted the 
origin from the core of the boojum by a distance a to the centre of the domain, i.e., the centre of the circle 
of radius ro and < < 2ir is the angle to the polar axis, 



Such a parametrization only includes shapes with no overhangs, i.e., the coordinate r(0) is a single-valued 
function of 9. Further the parametrization is smooth and so does not include shapes with cusps, as in [6, 25]. 
Since the area of the domain has to be A, we must have 



A perturbative analysis [171 ] involves perturbing the domain boundary about the circle and the texture about 
the optimal texture (virtual boojum, annuli or hedgehog). Such an analysis, carried out in [l7| for the virtual 
boojum, indicates that the circle is unstable to the n = 2 mode when a > (To, below a critical value of KjaR. 

We complement the linear stability analysis by a variational calculation. The reason we do this is because 
(a) linear stability analysis typically underestimates phase boundaries and (b) linear stability analysis does 
not give the (new) stable configuration in the regime of instability of the circle. To address these issues, we 
parametrise both the boundary and texture and determine the lowest energy configurations of the texture 
and the domain shape using a variational scheme. 

Again we set K = uq = 1 to fix the units of length and energy. We now minimise the energy as a function of 
variational parameters {TV, a, {c„}, {a„}} for a given value of R and a. The optimal value of TV is, as before, 
N = 2. We also find that including terms in the texture with coefficients {c„>i}, reduces the minimum value 
of E by at most 1%. So we parametrise the texture in the simpler form 



B. Noncircular domains 




(13) 




(14) 



<j> = 20 + c lV 



(15) 



where y is the distance of any point from the axis of the domain. As a is increased from 0, the core of the 
boojum approaches the centre of the domain but the domain remains almost circular. 
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FIG. 5: Stretching action of strong anisotropic line tension - domain conformation for (A) a = 6.5, (B) a — 9.5. (C) 
variation of 1/a (thin line) and the prolateness (/ — b)/b with a for a domain of effective radius R — 1. 



As cr is raised beyond a threshold a* , the core begins to recede from the centre and the domain bulges out 
at the equator (Fig. 5) and flattens near the poles. Indeed we find that the texture throughout the (R, a) 
plane is either a virtual boojum or a hedgehog defect. For large values of a, but still in the boojum phase, 
it is a good approximation to set <p — c±y and take the domain to be a rectangle bracketted at the equator 
by two semicircular caps (Fig. 5). 

Such rectangular domains have been observed [2(| on giant unilamellar vesicles (GUV's) consisting of a 
binary lipid mixture, a saturated component (e.g., DPPC) and an unsaturated component (e.g., DOPC), at 
temperatures below the gel transition of the saturated lipid (e.g., 42° C for DPPC). Following a temperature 
quench, domains of the saturated lipid, observed using fluorescence microscopy, were found to nucleate in 
the background of the unsaturated lipid and grow as rectangular domains. We believe that the tilt of DPPC 
in the gel phase gives rise to an anisotropic line tension, and can account for the shape of these domains, 
consistent with our analysis. 



C. Chiral domain 



We now turn to a description of the phases with nonzero bulk chirality k c . In this case, both the texture 
and the shape of the boundary may assume chiral shapes. For the moment however, let us assume that 
the boundary of the patch is circular. As before we set k\ = k% = 1 and o~q = 1. Further without loss of 
generality we take k c > 0. To highlight the effects of bulk chirality, we have set the anisotropic line tension 
(7 = 0. It will be clear that the effects of bulk chirality dominate the nature of the texture. 

Let us rewrite the Frank expression as 

E = 2ttR+ I - (divm + curlm) 2 + (k c - l)(divm)(curlm) (16) 
J A 2 

It is clear, and we have checked this explicitly, that for low enough bulk chirality, \k c \ < 1, the Uniform phase, 
with a circular domain and uniform m, is the lowest energy state. Indeed had we included the anisotropic 
line tension, we would have seen that the optimal achiral domains described in Fig. 4 would be the lowest 
energy configurations for small enough \k c \. 
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FIG. 6: (A) Pure divergence (B) Pure rotational (C) Spiral defect. 



As seen in (|16j) . increasing the magnitude of k c , makes the optimal texture more wound up : the texture 
prefers to have a very high curl and a divergence equal and opposite to the curl. The optimal texture is 
neither a pure divergence nor purely rotational, but an (Archimedes) spiral with an angle of opening a — tt/4 
with respect to the local radial (Fig. 6). 

In polar coordinates, with the origin at the centre of the circular domain, the spiral hedgehog is described by 
an m with constant radial and tangential components everywhere, m r = cos a and mg = sin a (m^+m^ = 1). 
Using 

divm= !^ (17) 
r 

CU rlm=^ (18) 
r 

we may compute the energy of this spiral configuration, 

1 R 

E = 2n(r c + R) - - (k c sin 2a - 1) In — + e c (19) 

2 r c 

where r c is the radius and e c is the energy of the defect core. The optimum value of a is 7r/4 when \k c \ > 1, 
making m r = rn$ — —l/y/2. The spiral hedgehog phase clearly gives a nonzero value for the chiral 

order parameter C (jlip . 

We see that the effect of the fc c -term is felt most strongly at the core of the spiral defect, falling off inversely 
as the square of the distance from the core. Can we find a texture whose chiral strength is large not just at 
one point but over the entire domain ? If so, this would surely be a candidate for the optimal texture. 

To answer this question, we have carried out a Monte Carlo simulation of 3055 particles carrying 0(2) spins 
on a triangular lattice with a hamiltonian obtained by discretising E. In order to search for the lowest energy 
configuration, we had to resort to simulated annealing from (ks T) -1 = 0.1 — * 300 starting from a variety 
of initial conditions. The energy functional is augmented by a higher order term 7 ((divm) 4 + (curlm) 4 ) 
to provide a cutoff to the spatial variations of m. Fig. 7 shows the optimal texture obtained for a = and 
k c = 12.75 — the texture is best described as a chiral tweed with a stripe size I* = 0.01. 

This texture may be easily parametrised variationally. First place the origin of coordinates on the boundary 
of the domain with the x-axis along a diameter, and parametrise the texture by 

x x 
m x — I cos y| , m y — -| siny | (20) 

The texture consists of stripes of width irl parallel to the y-axis. This parametrisation is shown in Fig. 7 for 
comparison with the optimal texture obtained in the simulation. This texture confers a net chiral strength 
to the domain — (V ■ m)(V x m) has the same sign everywhere, though the sign of individual terms, V • m 
and V x m, varies from one stripe to the next. The energy of the chiral tweed domain is 

E = 2ttR-^J dx^j2xR-x 2 \An^\ + ^ r (21) 



FIG. 7: (a) Close up of the chiral tweed texture generated by a monte carlo simulation, and (b) its continuum 
representation by a mathematical formula. On the shaded (clear) portions div m is positive (negative) and curl m is 
negative (positive) . This shows the texture within a unit cell ; the pattern repeats periodically to form the Tweed 
phase. 




FIG. 8: Phases of the tilt texture domain with bulk chirality and a = : (1) uniform phase (2) spiral defect (where 
we have taken e c = and r c = 0.005) (3) chiral tweed (stripe width I* = 0.01). 



The variational parameter I approaches zero to minimise this energy — higher order derivative terms in the 
hamiltonian, such as 7((divm) 4 + (curlm) 4 ), would however cutoff this monotonic decrease at some scale 
I*. Note that the energy decreases rapidly as k c increases (Fig. 9). 

The phase diagram of this chiral domain is given in Fig. 8, together with a demonstration of discontinuous 
transitions as seen by the crossing of the energy branches, Fig. 9. 




FIG. 9: Variation of the energy branches of the (1) Uniform (2) Spiral Hedgehog and (3) Tweed phases as a function 
of k c at R = 1 indicating discontinuous transitions. Notice that the energy decreases rapidly in phase (3) as k c 
increases. 
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FIG. 10: Energy E^ n \ of a domain split into n equal parts, depends on its total area, A = nR 2 . 
n = 1 and 2 branches : the domain prefers to split in two for a range of R values (aor c = 1.5, e c = 



Here we show the 
= 5, k c = 17). 



IV. DOMAIN SPLITTING 



So far our restriction to the strong segregation limit has also assumed that the planar domain of area A 
does not break up. If the chiral strength is large enough, k c 3> 1, the texture might prefer to maximise the 
number of spiral defect points. This could induce domain splitting. To study the conditions under which 
such breakup is favourable, we calculate the energy of n circular domains of equal area, each bearing the 
same spiral texture and compare it to a single circular domain with the same total area and texture. The 
total energy of this configuration is 

£jW = 27rcr \ALR - nir (k c - 1) In ^H^L + n€c (22) 



For small values of A, a single domain E^> has the least energy. As A increases, becomes smaller 
than E^ : chirality in the bulk wins over interfacial energy (Fig. 10). As A increases further, multidomain 
splitting is favoured. At even higher values of A, however, domain splitting is prohibited by a large interfacial 
energy cost; thus there is a window of areas for which chirality induces domain split up. This tendency to 
split holds when k c is large enough; for a fixed value of crrjTc, there is a critical k c beyond which chirality 
induced splitting would manifest. 

This chirality induced splitting has implications in the nucleation and growth of tilt-ordered domains in 
a fluid substrate, as in binary lipid mixtures where one of the components has a tilt and is chiral or in 
the nucleation of Sm-C* domains in a Sm-A film (since the SmA-SmC* transition can be first order [271]). 
For instance, consider a GUV composed of a binary lipid mixture of a saturated (minority component) 
and unsaturated lipid species, quenched below the gel transition of the saturated lipid. Domains of the gel 
phase spontaneously nucleate in the liquid-disordered phase. Growth initially occurs via coarsening; as the 
domains get larger, they undergo random brownian diffusion on the surface of the GUV, coalescing when 
two domains encounter each other. If however the molecules in the gel phase have a tilt and chirality, then 
coalescence could be prevented by the above mechanism when the domains reach a size R*/y/2 (Fig. 10). 
Further growth of the domain would be halted by this chirality induced domain repulsion. A similar feature 
should be observed in the nucleation and growth of Sm-C* domains in a Sm-A film. 
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In this section we estimate the various parameters entering the domain free-energy for typical systems under 
study. We then compare the Frank energy scale to dipolar energies and provide some justification for why 
we have ignored the latter. 



The values of the Frank constants entering the domain energy functional ©, may be obtained from the 
corresponding values in bulk liquid crystals. In a cholesteric liquid crystal, the director field ft, describing 
the locally averaged molecular (long)-axis, describes a helical conformation about a fixed ordering axis. The 
energy density of this conformation is given by the Frank expression [l5j 

£ = (div n) 2 + ^ (ft • curl ft + q) 2 + (n x curl nf (23) 
where 2ir/q is the pitch of the helical conformation, much larger than the molecular size. 



Comparing the energy density of the planar tilt-domain ^ with that of a cholesteric above, suggests the 
following correspondence : k\ ~ aK\ , and /C2 ~ a K% , where the length scale a is of the order of the thickness of 
the surface bearing the domain. In addition, £ has a term linear in curl ft, suggesting that k c ~ aK2- Knowing 
the values of the Frank constants for the bulk cholesteric, K\ ~ K2 ~ K3 w 10 -11 N, and taking a w 10A, 
we arrive at an estimate for the Frank coefficients in the energy functional (J3J), k\ ~ k^ ~ k c as 10~ 20 J. 
The magnitude of the chiral coefficient k c could be enhanced by the addition of chiral impurities. Note that 
these values are almost an order of magnitude larger than fc^T at 30° C, and so we may ignore the effect of 
thermal fluctuations on the texture. 



Estimates of the line tensions may be obtained from studies of shape and texture changes of tilt domains 
in Langmuir monolayers using Brewster angle microscopy [5| and shapes of lipid domains in GUV's using 
fluorescence microscopy [9(. These methods give an estimate of Co ~ 4 x 10~ 13 N for the isotropic line 
tension, and a value of 3.6 x 10 4 m _1 for the ratio a/K. With our estimate of 10 -20 J for the Frank constant 
K, we arrive at an anisotropic line tension a as 3.6 x 10^ 16 N. 



With these parameter estimates, we now argue post facto that long-ranged dipolar interactions are signifi- 
cantly weaker than the Frank contributions. Typical phospholipids have a dipole moment \p\ « 1 debye Q 
directed roughly parallel to the plane of the membrane. The strength of dipolar interactions of neighbour- 
ing tilted phospholipids is of the order \p\ 2 /a 3 where a w lOA is the separation of the lipid dipoles. This 
energy is of the order of 10~ 22 J, as a result, dipolar interactions can not challenge the order imposed on 
m by the Frank energy. However if the molecules are charged, then the strength of electrical interaction 
may be considerably enhanced \p\q/a 2 w 10~ 20 J for q = 1 coulomb and a as 8A, giving rise to significant 
chargc-dipolc or dipolc-dipolc interactions. Colloidal systems with substantial dipolar effects are known to 
exhibit string-like aggregates [2lj . 



VI. DISCUSSION 



In this paper we have discussed the equilibrium shapes and textures of a single domain consisting of molecules 
described by a tilt embedded on a flat substrate containing non-tilt molecules. In addition to having a tilt, 
the molecules constituting the domain may be chiral. The interplay between tilt, geometry and chirality 
give rise to a rich variety of shapes and textures separated by discontinous transitions. We have arrived 
at our results by a combination of Monte Carlo simulations and exact variational calculations. In addition 
our variational ansatz allows for multiply connected domains. Our key results : (i) In the achiral limit, 
optimum textures include virtual boojum, annular and hedgehog phases. A novel feature is the occurence of 
spontaneous chiral symmetry breaking to a chiral annular phase, (ii) When the domain shapes are allowed 



to deviate from circularity, then the tilt induced anisotropic line tension typically give rise to domain shapfes 
which are flattened and elongated, (iii) Chirality produces spiral defects and an unusual chiral tweed phase, 
(iv) Chirality induces large enough domains to breakup into smaller domains leading to a limiting domain 
size (keeping all other parameters fixed) and preventing coalescence. 

On symmetry grounds, one may argue [lil . [l5T ] that the tilt phase of chiral molecules having a permanent 
dipole moment is ferroelectric. Thus an externally applied electric field E can couple to the polarization 
vector P, via an E ■ P coupling. The polarization vector P is perpendicular to both the tilt m and the normal 
to the plane z |l5j, i.e., P = az x rh = a(s'm </>, cos <f>, 0). The total energy ©, (0J must be augmented by 
—a(E x sin <fi + E y cos0), in addition to an electrostatic self energy proportional to P 2 . In the low chirality 
limit k c < 1, this field imposes a preferred orientation for m; by changing the direction of the field, one may 
switch between the two textures [16[ . In addition, strongly chiral systems admit a pseudo-scalar coupling of 
the form (m • Vm • P), which is simply (V x m) z [15]; this leads to a switch in the handedness of the spiral 
texture on changing the direction of the external field. 

While the results obtained here are applicable to a variety of systems featuring tilt domains of finite size 
embedded in a 2d substrate, our main interest is in the shapes and dynamics of domains in a multicomponent 
lipid mixture. Giant unilamellar vesicles composed of saturated and unsaturated lipid components, typically 
exhibit a wide coexistence regime between a liquid disordered phase (enriched in the unsaturated component) 
and a gel phase (enriched in the saturated lipid). Following a quench into this coexistence region, domains 
of the minority component nucleate and grow in the majority phase. If the saturated lipid is chiral and has 
a tilt in the gel phase, then we believe that the analysis carried out here has some relevance to the shapes 
and dynamics of nucleating domains [2(| H^] . 
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